Controllable spin-dependent transport in armchair graphene nanoribbon structures 
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Using the non-equilibrium Green's functions formalism in a tight binding model, the spin- 
dependent transport in armchair graphene nanoribbon (GNR) structures controlled by a ferro- 
magnetic gate is investigated. Beyond the oscillatory behavior of conductance and spin polarization 
with respect to the barrier height, which can be tuned by the gate voltage, we especially analyze the 
effect of width-dependent band gap and the nature of contacts. The oscillation of spin polarization 
in the GNRs with a large band gap is strong in comparison with 2D-graphene structures. Very 
high spin polarization (close to 100%) is observed in normal-conductor/graphene/normal-conductor 
junctions. Moreover, we find that the difference of electronic structure between normal conductor 
and graphene generates confined states in the device which have a strong infiuence on the transport 
quantities. It suggests that the device should be carefully designed to obtain high controllability of 
spin current. 

PACS numbers: 75.75.+a, 72.25.-b, 05.60.Gg, 73.43.Jn 



I. INTRODUCTION 



Graphene, a monolayer of carbon atoms packed into a 
two dimensional (2D) honeycomb lattice, has attached a 
great amount of attention from both experimental and 



since it was isolated and 
[^. It is a basic building 



theoretical points of view 
demonstrated to be stable 
block for graphite materials of all other dimensionalities, 
e.g., it can be wrapped up into OD fuUerenes, rolled into 
ID nanotubes, or stacked into 3D graphite. Due to its 
unique electronic properties, i.e., its conduction electrons 
behave as massless Dirac fermions [1, 3 j a lot of inter- 
esting phenomena such as the finite conductance at zero 
concentration [3J , the unusual half integer quantum Hall 
effect Q , and the Klein tunneling Q have been observed 
in the graphene and theoretically discussed in the frame- 
work of the massless fermion Dirac's model 3, 6]. 

Finite width graphene strips, which are referred as 
graphene nanoribbons, have been also studied attentively 
0, B S E [HI- It has been shown that GNRs of var- 
ious widths can be obtained from graphene monolayers 
using patterning techniques . The transport properties 
of the perfect GNRs are expected to depend strongly on 
whether they have zigzag or armchair edges Q. In the 
framework of the nearest neighbor tight binding (NNTB) 
model, the GNRs with zigzag edges are always metal- 
lic while the armchair structures are either semiconduct- 
ing or metallic depending on their width. In the zigzag 
GNRs, the bands are partially flat around the Fermi en- 
ergy [E = Q eV], which means that the group velocity 
of conduction electrons is close to zero. Their transport 
properties are dominated by edge states . In the GNRs 
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with armchair edges, the bands exhibit a finite energy 
gap in the semiconducting structures or are gapless in 
the metallic ones |8j. However, ab initio studies have 
demonstrated that there are no truly metallic armchair 
GNRs (see in Ref. ^ and reference therein). Even 
for the structures predicted to be metallic by the NNTB 
model, a small energy gap opens, thus modifying their 
behavior from metallic to semiconducting. In general, 
the group velocity of conduction electrons in armchair 
GNRs is high, e. g., it is constant and equal to about 
10^ m/s |3] in those metallic structures. The transport 
properties of ideal zigzag and armchair GNRs are thus 
very different. In the current work, we pay our attention 
only to the ribbons with armchair edges. 

Experimentally, electronic transport measurements 
through a graphene sheet usually require contacts to 
metal electrodes, e.g., see an illustration in Ref. p^ . 
When tunneling from metal reservoir to graphene oc- 
curs in a large area, the contact becomes Ohmic 
and the area under the contact forms a substance 
which is a hybrid between graphene and normal 
metal 13;]. Depending on the nature of this sub- 
stance, the system can be appropriately considered 
as a graphene/graphene/graphene (GGG) structure or 
a normal-conductor/graphcne/normal-conductor (NGN) 
junction whose contacts can be modeled by honeycomb 
or square lattices, respectively. The ballistic transport 
through the NGN systems has been investigated system- 
atically in Refs. [IlQIIl. 

Beside their interesting electronic transport properties, 
due to very weak spin orbit interaction 16], which leads 
to a long spin flip length Ijim) [17| . the graphene 
- based structures also offer a good potential for spin- 
polarized electronics. Actually, graphene is not a natural 
ferromagnet. However, recent works have shown that fer- 
romagnetism and spin-polarized states can be introduced 
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in graphene, e. g., by doping and by defect [l9|, |20|, |2l| 
or by applying an external electric field (2^ . Especially, 
Ifaugen et. al. [13] suggested that the ferromagnetic 
correlations can be created in graphene by the so-called 
proximity effect. The exchange splitting induced by de- 
positing a ferromagnetic insulator EuO on the graphene 
sheet was then roughly estimated to be about 5 meV. 
The effect has been also experimentally demonstrated 
[13, [ill- Motivated by these features, some other works 
have predicted and discussed the controllability of spin 
current by a ferromagnetic gate in 2D-graphene struc- 
tures [U |2^ [2^ . The spin current was found to be an 
oscillatory function of the potential barrier, which can 
be tuned by the gate voltage, and its amplitude is never 
damped by the increase of the width and the height of 
barrier [24| . However, the spin polarization is not so 
high, e. g., its maximum value is just about 30%. In ad- 
dition, the other spin-dependent properties of graphene 



spin Hall effects 
331 have been also 



such as spin field effect transistor 127 
[H [13, spin valve effects [H, M, m 
investigated extensively. Especially, the giant magneto- 
resistance has been explored theoretically and discussed 
in the structures of armchair 3l| and zigzag [33 | GNRs 
connecting two conducting electrodes. The authors pre- 
dicted that it can reach the high value of 100% [s^l • 

In this article, we are interested in the possibilities of 
electrically tunable of spin current in single ferromag- 
netic gate armchair GNR structures. By investigating 
the physics of spin polarized transport in these struc- 
tures, we would like to derive some simple scaling rules 
in order to tend to a high tunable spin polarized cur- 
rent. The study is focussed on the role of the ribbon's 
energy band gap and the different types of leads (either 
graphitic or normal-conducting). In the NGN systems, 
since the strength of the device-to-contact coupling and 
the device length are important parameters, their influ- 
ence on the control of spin current is also carefully in- 
vestigated. The paper is organized as follows. Section 
2 is devoted to the description of the model and main 
formulas based on the non-equilibrium Green's functions 
formalism (NEGF). In section 3, the numerical results 
are presented and discussed. Finally, a brief summary is 
given in section 4. 



II. MODEL AND FORMULATION 

The considered structures consist of an armchair GNR 
coupled with two electrodes which may be described 
either by graphitic (Fig. 1(a)) or normal-conducting 
(Fig. 1(b)) leads. In the simplest consideration, the 
normal-conducting leads are modeled by square lattices 
[13I [1^ [lH . A ferromagnetic gate is assumed to create 
a potential barrier which controls the Fermi level locally 
and to induce an exchange splitting into the device. To 
model the structures, we use the single band tight bind- 
ing Hamiltonian 
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FIG. 1: (color online) Schematic illustration of the consid- 
ered armchair GNR structures with the number M of car- 
bon chains between two edges: (a) graphitic and (b) normal- 
conducting leads. The latter ones are modeled by square lat- 
tices. A magnetic gated insulator is deposited to create a 
spin-dependent potential barrier in the center of device. 



where Hl,fi are the Hamiltonian of the left and right 
leads, respectively; Hd is the Hamiltonian of the device; 
He describes the coupling of the device to the leads. The 
Hamiltonian terms in Eq. (1) can be written as 



Ha 

Hd 
He 



-tc (cL,^aid,^ + /i-C-) 

a={L,R} {ia-jd),<y 



where the operators cj^ ^ {ci^^a) and ^ {o,id,<y) create 
(annihilate) an electron with spin a in the electrode a 
and the device region, respectively. The sum over carbon 
atoms (i, j) is restricted to the nearest neighbor atoms, i, 
and tc stand for the hopping parameters in the device, 
the lead and at the coupling interface, respectively. 
is the on-site energy of the leads which acts as a shift in 
energy. The device spin-dependent on-site energy £^^,0- is 
modulated by the gate voltage 



Ua-ah 




in gated region 
otherwise 



(3) 



H = Ht 



Hr 



H„ + Hr 



(1) 



Here, Ug denotes the potential barrier height, h is the 
exchange splitting and cr = ±1 describes the up/down 
spin states. 

Since no spin flip process is considered here, (1) can 
be decoupled into two linear spin-dependent Hamiltoni- 
ans Ha^ and the transport is easily considered using the 
NEGF formalism. For each spin channel cr, the retarded 
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Green's function is defined as 



Gl (E) 



E + iO+ - HD,a - Sl - 



(4) 



where, describes the retarded self energy matrices 
which contain the information on the electronic struc- 
ture of the leads and their coupling to the device. It can 
be expressed as = TD,a9aT^a,D where f is the hopping 
matrix that couples the device to the leads, ga are the 
surface Green's functions of the uncoupled leads, i.e., the 
left or right semi- infinite electrodes. The surface Green's 
functions and the device Green's functions are calculated 
using the fast iterative scheme [ssj and the recursive al- 
gorithm [U, respectively. 

Within the model described by Eq. (1), the transport 
is considered to be ballistic and the conductance through 
the device is calculated using the Landauer formalism 
[35| . The spin-dependent conductances Ga at the Fermi 
energy Ep are obtained from the transmission function 
Ta{E), such that 



and 



Ga (Ef) = -T, {Ef) 



Ta (E) = Tr FiG^r^G^ 



(5) 



(6) 



Here, G[^ (= G^J) denotes the advanced Green's function. 
The tunneling rate matrix rL(_R) for the left (right) lead 
is obtained from 



(7) 



where S° (= T,^) is the advanced self energy. Finally, 
the spin polarization is determined by 



P = 



G] + Gi 



(8) 



In addition, the local density of states (LDOS) at site j 
can be also directly extracted from the retarded Green's 
function as: 



LDOS(,) = — ImG'^(j,i) 



(9) 



By using the recursive algorithm described in Ref . [3J] , 
the size of matrices equals the number M of carbon 
chains between two edges. Thus, the cost of calculations 
is only linearly dependent on the device length. 



III. RESULTS AND DISCUSSION 

Using the above formalism, we investigate the 
spin-dependent transport in the considered structures. 
Throughout the work, we set t = 2.66 eV [s^ in the 
graphitic regions and assume that is equal to t and 
tc < t in the NGN junctions [H, Our consideration 
is restricted to the low-energy regime when E <^t. 
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FIG. 2: (color online) (a) Oscillation of conductance vs 
the barrier height Ua in the GGG structures with different 
widths: M = 21 (dashed) and 23 (solid line), (b) illus- 
trates the transmission coefficient calculated from eq. (10) 
for different modes 9j. Other parameters are: Lq ~ 42.5 nm, 
Ef = 300 meV and /i = meV. 



A. Gate control of spin current in GGG structures 

As mentioned above, the gate voltage creates a poten- 
tial barrier in the device. In the armchair GNRs, the elec- 
tronic properties of the structures with AI ^ 3n + 2 and 
M = 3n -|- 2 (n is an integer) are significantly different, i. 
e., there is a finite energy band gap in the former struc- 
tures while it is negligible in the latter ones [TO] • Respec- 
tively, we display in Fig. 2(a) the conductance as a func- 
tion of the barrier height Ug for the structures: M = 21 
(dashed) and 23 (solid line). Here, the gated region is as- 
sumed to be nonmagnetic, i. e., /i = meV. First of all, 
the results show that the conductance has an oscillatory 
behavior with respect to Ug- This phenomenon has been 
observed in the 2Z)-graphene structures and explained as 
a consequence of the well-known Klein's tunneling [6l.l23|. 
In the framework of the Dirac's description, the conduc- 
tance peaks have been demonstrated to be essentially 
due to the resonance or the good matching of electron 
states and confined hole states outside/inside the bar- 
rier region [2^, respectively. Those states, in this model, 
correspond to electron states in positive and negative en- 
ergy bands. Note that due to the finite ribbon width, 
the transverse momentum is quantized into a set of dis- 
crete values. Practically, the oscillation of conductance 
can be seen clearly from the analytical expression of the 
transmission coefficient for a given transverse momentum 
mode fc^ (see the calculation in Ref. [131 )■ In the limit 
of low energy, it can be rewritten as 



T = 



cos2 0, sin^ , 



cos^ 9j sin^ (j) + (sin 9j + cos 0) sin^ (k^LG) 



(10) 



where 9^ 



tan ^{kl/k^ 



and 



[{t-YFkl)s\n{iakl/2)] / [Ug- E) with v_f = 3a</2 
and a is the C — G bond length, kx (fc^) denotes 
the longitudinal (transverse) momentum, which is the 
deviation of the momentum k from the zero energy 
point, outside the barrier and k\ is the longitudinal 
momentum inside the barrier. The energy dispersions 
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outside/inside the barrier are, respectively 

E - VFv/fc^+fcf (11) 
E-Ug = ypkl) sin^ ^ + v|/fcf (12) 

Accordingly, the transmission and then the conductance 
have their maximum (or minimum) values when k^Lc is 
equal to rmr (or (m + l/2)7r) for any integer m. In the 
limit oi E Uq ^ i, the eq. (12) can be rewritten as 
Ug — E K, Ypk!^ and the period of oscillation is defined 
by Up = yp'kILg-, which coincides with that in Refs. 
[2J, [2^. In general cases, when the relation between 
Vg and k is nonlinear, the period can be approximately 
expressed B&Up — YqIt/Lg with v^ < vf- For instance, 
Vg « 0.74vi^- and 0.89vi? for M = 21 and 23 presented in 
Fig. 2(a), respectively. 

As a consequence of different electronic structures, the 
Fig. 2(a) also shows that the conductance for the case 
of M = 21 (large energy gap) oscillates strongly in com- 
parison with the other (M — 23). This can be easily 
understood by considering the behavior of transmission 
coefficient for different energy gaps. As seen in Eq. (11), 
the energy gap is enlarged when increasing fc^ (or 9j). 
So that, from the Eq. (10) and the Fig. 2(b), we see 
that when the energy gap is larger (increasing 9j), the 
oscillation of transmission is stronger. It leads to the 
different behaviors of conductance shown in Fig. 2(a). 
Similarly, it is shown that the oscillation of conductance 
in the structure with M = 21 is also stronger than that 
in the 2Z?-graphene structures (see in Fig. 2(a) of Rcf. 
(23|). where the gap is truly zero. 

In eqs. (10-11), the mode 6'o = (or ky = 0) corre- 
sponds essentially to the normal incident mode whose en- 
ergy dispersion is gapless and linear in the 2_D-graphene 
structures. When considering the behavior of transmis- 
sion coefficient, we find an important feature: it is not 
uniformly equal to unity, but a function of the barrier 
height even in the case of 6* = ^o- This differs from the 
prediction of Klein's paradox obtained by using Dirac's 
description in graphene structures Q as previously dis- 
cussed in Ref. [38]. Practically, the transmission coeffi- 
cient (eq. (10)) approaches the simphfied expression (4) 
in Ref. [a| only in the limit oi E <^Ug <.t- 

Now we investigate the behavior of spin polarization in 
the ferromagnetic gate structures. The exchange split- 
ting h is chosen to be 10 meV, which can be achieved 
experimentally [l3, [l3|- Since no spin flip process is 
considered, the exchange splitting just shifts the conduc- 
tance of each spin channel relatively to the other. The 
spin polarization therefore behaves as an oscillatory func- 
tion of Ug as shown in Fig. 3. Similar phenomena in 
the 2Z?-graphene structures have been also observed and 
discussed in Rcf. 0, [2^, |2^. It was shown that the 
oscillation of spin current is never damped with the in- 
crease of the width and the height of barrier and the 
spin polarization can be reversed by changing the gate 
voltage. Actually, the amplitude of P depends primarily 
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FIG. 3: (color online) (a) Spin polarization as a function 
of the barrier height Ug for the same structures as in Fig. 
2(a). (b) shows an example of the effects of the different rib- 
bon widths on the spin polarization: AI = 21 (dashed), 27 
(dashed-dotted) and 33 (solid line). Everywhere Lg = 42.5 
nm, Ef = 300 meV and h = 10 meV. 

on the phase coherence/decoherence of the oscillation of 
spin - dependent conductances, i. e., it has the maxi- 
mum/minimum value when the gate length (or the bar- 
rier width) Lg is equal to a half-integer/integer of Lh 
with Lh ~ Vg7r/2/i, respectively. Hence, the gate con- 
trol of spin current can be modulated by changing Lg, 
i. e., it leads to the beating behavior of P similar to 
that shown in Fig. 5(c) of Ref. [2^. Furthermore, as a 
consequence of the behavior of conductance presented in 
Fig. 2(a), Fig. 3(a) also shows that the oscillation of P 
in the GNRs with a large energy gap is very strong in 
comparison with the others. For instance, the amplitude 
of P is about 65% for M = 21 while it is only few per- 
cents for AI — 23 or has a maximum value of 30% in the 
2£'-graphene structures [l^, H^H^. However, since for 
M 7^ 3n -f 2 the energy gap decreases when increasing 
the ribbon width, the oscillation of conductance and P 
in those structures is gradually weaker. The transport 
quantities for the GNRs in the limit of infinite width ap- 
proaches those for the 2Z3-graphene structures, where the 
continuum Dirac's description is valid. To illustrate this 
point, we display in Fig. 3(b) an example of the effect 
of different ribbon widths on the spin polarization in the 
structures with M ^ 3n-|-2. Indeed, when increasing the 
ribbon width, the amplitude of P decreases and becomes 
closer to that in 2D-graphene, i. e., it is only about 35% 
for M = 33. 



B. Effects of normal-conducting leads 

In this section, we consider the spin transport in NGN 
junctions. First, we focus on the possibilities of obtaining 
high tunable spin current when replacing the graphitic 
leads by the normal-conducting ones. Second, we an- 
alyze the sensitivity of transport quantities to different 
parameters as the device length and the Fermi energy. 

In Ref. [l3|, Schomerus compared the resistances of 
NGN junctions and GGG structures, and found the dual- 
ity between graphitic and normal-conduncting contacts. 
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FIG. 4: (color online) Comparison of conductance (a) and 
LDOS (b) in different structures: graphitic (dashed) and 
normal-conducting (solid lines, tc ~ t) leads. Everywhere 
M = 21, Ld = 51 nm, La = 42.5 nm, Ef = 300 meV and 
ft = meV. 



He has shown that identical transport properties arise 
when the graphitic leads are replaced by quantum wires 
and the difference between the results obtained in those 
structures is only quantitative. On this basis, we plot 
in Fig. 4(a) the conductance as a function of Uq in the 
NGN junction in comparison with the GGG structure 
for the case of Ljj = 51 nm and Ep — 300 meV. Since 
we found that the results depend weakly on the hopping 
energy in the leads, ti is chosen to be equal to t for sim- 
plicity. Qualitatively, Fig. 4(a) shows that the oscillation 
of conductance seem to be unchanged in its phase and pe- 
riod when changing the leads. Quantitatively, the oscil- 
lation in the NGN junction is stronger than in the GGG 
structure. This can be explained clearly by the effect of 
replacing the graphitic leads by the normal conducting 
ones on the picture of bound states in the barrier region. 
These states, in the framework of the Dirac's description, 
have been considered as confined hole states in the 2D- 
graphene structures [l^l. In Fig. 4(b), we display the 
LDOS (see the right axis for graphitic and the left axis 
for normal-conducting leads) at the first site of barrier 
region with respect to Ug- The results are understood 
that the peaks of LDOS occur when the Fermi energy 
corresponds to any bound state. Obviously, the oscil- 
lation of LDOS (or the quantization of bound states in 
the barrier) in the NGN junction appears stronger (with 
higher peaks) than that in the GGG structure. It is the 
essential origin of the different behaviors of conductance 
as shown in Fig. 4(a). 

Now we turn to the behavior of spin current in ferro- 
magnetic gate NGN junctions, i. e., /i = 10 meV. In Figs. 
5(a) and 5(b), we display the comparison of spin polar- 
ization in the NGN junctions and the GGG structures. 
Due to the different behaviors of conductance shown in 
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FIG. 5: (color online) (a,b) Comparison of spin polarization 
in the different structures: graphitic (dashed) and normal- 
conducting (solid lines, tc ~ t) leads. (c,d) The spin po- 
larization in the latter one with different coupling strengths; 
tc ~ t (dashed), 0.8t (dashed-dotted) and 0.6t (solid lines). 
The ribbon widths are M = 21 (a,c) and 23 (b,d). Other 
parameters are Ld = 51 nm, Lq = 42.5 nm, Ep = 300 meV 
and h = 10 meV. 



Fig. 4(a), the amplitude of P in former structures is re- 
markably larger than that in latter ones. Particularly, 
when changing the leads, it increases from 65% to 81% 
(see in Fig. 5(a)) and from 1% to 50% (see in Fig. 5(b)) 
for M — 21 and 23, respectively. Moreover, in the NGN 
junctions, the possibility of obtaining high tunable spin 
current is more impressive with decreasing the strength 
of the dcvice-to-contact coupling, which is characterized 
by the hopping energy tc- In Figs. 5(c) and 5(d), we 
plot the obtained results for three cases: tc = t (dotted), 
0.8i (dashed) and 0.6t (solid lines). Actually, the trans- 
port in these structures depends strongly on the prop- 
erties of junctions and therefore ow tc- A smaller tc 
corresponds to a higher contact resistance We find 
that the quantization of bound states in the barrier re- 
gion is stronger when decreasing tc (not shown), which 
leads to a stronger oscillation of the transport quanti- 
ties with respect to the barrier height. Indeed, even in 
the case of M = 23, the spin polarization can reach a 
very high value of 86% for tc — 0.6t (see in Fig. 5(d)). 
More impressively, in the cases of A/ = 21, it can tend to 
100% by reducing tc (see in Fig. 5(c)). A similar feature 
(giant magneto resistance) has been also predicted in the 
structures of GNRs connecting two conducting electrodes 

MM- 

Practically, the features discussed above depend 
strongly on the parameters of the NGN junctions, such 
as the device length and/or the Fermi energy. It results 
from the fact that the charge transport can be confined 
in the device by two normal-conductor/graphene junc- 
tions. It leads to an additional resonant condition con- 
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FIG. 6: (color online) (a) LDOS illustrating the existence 
of confined states in the device and (b) conductance in the 
NGN junctions as a function of the Fermi energy for different 
device lengths: Ld ~ 51 nm (solid) and 102 nm (dashed 
Unes). (c) shows the dependence of energy spacing of the 
confined states on the inverse of device length. Everywhere 
M = 21, Lg = 42.5 nm, tc = 0.8t and Ua = h = meV. 

trolling the transport picture beside the transmission via 
the bound states in the barrier. Indeed, the existence 
of such confined states is demonstrated clearly from the 
behavior of LDOS shown in Fig. 6(a). In this figure, to 
cancel the effects of bound states in the barrier, the gate 
voltage is not applied to the device. The energy spacing 
Es is estimated to be about 25.5 meV for Ld = 51 nm 
and 12.6 meV for 102 nm. It means that Es seems to 
be inversely proportional to the device length, i. e., as 
illustrated in Fig. 6(c). This implies an unusual quanti- 
zation of charges in the graphene-based structures, which 
is essentially different from the case of normal semicon- 
dutors wherein Es oc 1/Lj~i as previously discussed in 
Refs. [25;, |3S]. Due to such confinement, the transport 
quantities, such as the conductance and the spin current 
(not shown), have an oscillatory behavior also with re- 
spect to the Fermi energy in the considered region (see 
in Fig. 6(b)). Therefore, when a gate voltage is applied, 
there is a coexistence of bound states in the barrier and 
confined states in the device. They together respond for 
the resonant transport conditions of the structure. 

On the other hand, the gate controllability of spin cur- 
rent is principally due to the picture of bound states in 
the barrier (or Klein's tunneling) [2^1 • It arises a ques- 
tion about how the confined states in the device affect 
that picture. As shown in Figs. 4 and 5, the replacement 
of the graphitic leads (infinite Ld) hy normal-conducting 
ones (finite Ld) does not affect the period, but the ampli- 
tude of the oscillation. It suggests that, in the NGN junc- 
tions, the oscillation of conductance and spin polarization 
can be modulated in its amplitude while its period is un- 
changed when changing the device length. To examine 
this statement, we plot the conductance in Fig. 7(a) and 
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FIG. 7: (color online) (a) Conductance and (b) spin polar- 
ization P in the NGN structures as functions of the barrier 
height Ug for different device lengths: 68 nm (solid), 74 nm 
(dashed) and 79 nm (dashed-dotted lines). The oscillation of 
P versus the device length for different values ol Ug (c): 630 
meV (solid), 655 meV (dashed-dotted) and 681 meV (dashed 
line); and Ef (d): 250 meV (dashed-dotted), 300 meV (solid) 
and 350 meV (dashed line). Other parameters are AI — 21, 
Lg = 42.5 nm, Ef = 300 meV (in (a,b,c)), tc = 0.8t, 
Ug = 630 meV (in (d)) and ^ = 10 meV. 



the spin polarization in Fig. 7(b) as functions of the bar- 
rier height for different device lengths. From Fig. 7(a), 
we see that while the period is determined only by the 
gate length, the oscillation of conductance is modified by 
changing Ld, i. e., it is strong/weak when Ld = 68 (and 
79) or 74 nm, respectively. Consequently, the amplitude 
of spin polarization is dependent on Ld, i. e., it is about 
95% loT Ld^ 68 nm, 15% for 74 nm and 86% for 79 nm 
(see Fig. 7(b)). This feature is exhibited more clearly 
in Fig. 7(c) by three curves of spin polarization versus 
Ld for different barrier heights: Uq = 630 (dashed), 655 
(dashed-dotted), and 681.5 (soHd fine) meV. We see that 
the spin current has an oscillatory behavior and is sup- 
pressed completely at certain values of Ld- Obviously, 
this demonstrates that the amplitude of spin polarization 
exhibited in Fig. 7(b) is also an oscillatory function of 
the device length. Moreover, its period seems to be in- 
versely proportional to the Fermi energy, i. e., it is about 
27.8 nm for Ep = 250 meV, 18.3 nm for 300 meV and 
13.4 nm for 350 meV (see Fig. 7(d)). It is nothing, but 
a consequence of the resonant transport due to the con- 
fined states in the device. Hence, the gate control of spin 
current in the NGN junction can be modulated not only 
by the gate length Lq (see in the section A) but also by 
the device length Ld and/or the Fermi energy Ep. This 
implies that the structure should be carefully designed to 
obtain high controllability of spin current. 
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IV. CONCLUSIONS 

Using the NEGF method for quantum transport sim- 
ulation within a tight binding hamiltonian, we have con- 
sidered the spin-dependent transport in single ferromag- 
netic gate armchair GNR structures. The leads are mod- 
eled as either graphitic or normal-conducting. 

In the case of graphitic leads, it is shown that the con- 
ductance and the spin current behave as oscillatory func- 
tions of barrier height which can be tuned by the gate 
voltage. The oscillation of spin polarization in the rib- 
bon structures with a large energy band gap is strong 
in comparison with the 2£)-graphene structures. Es- 
pecially, the study has demonstrated that a very high 
spin polarization can be observed in the NGN junctions. 
It results from the fact that the quantization of bound 
states in the barrier (gated) region can appear very strong 
when using the normal-conducting leads. In this struc- 
ture, it is shown that the spin polarization increases and 
can tend to 100% with the decrease of the strength of 



the device-to-contact coupling. Moreover, we have also 
found the existence of confined states in the device by 
normal-conductor/graphene junctions. This confinement 
responds for an additional resonant condition beside the 
transmission via the bound states in the barrier. There- 
fore, the gate control of spin current in the NGN junc- 
tions can be modulated not only by the gate length but 
also by the device length and/or the Fermi energy. 

Our predictions may be helpful for designing efficient 
spintronics devices based on perfect armchair GNRs. 
However, some disorder effects, e. g. due to edge rough- 
ness, have been observed experimentally [3 and demon- 
strated to affect the transport properties [ll| of the GNR 
structures. Further work is needed to assess their influ- 
ence on the spin polarized properties discussed in this 
article. 
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